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Abstract. We study the single spin asymmetry (SSA) for the D-meson production A® in 
the pp collision, p^p — » DX, in the framework of the collinear factorization. Since the charm 
quark is mainly produced through the cc-pair creation from the gluon-fusion process, this is an 
ideal process to probe the twist-3 triple-gluon correlation functions in the polarized nucleon. 
We derive the corresponding cross section formula for the contribution of the triple-gluon 
correlation function to A® in p^p — > DX, applying the method developed for ep^ — > eDX 
in our previous study. As in the case of ep^ — > eDX, our result differs from a previous study 
in the literature. We will also present a simple estimate of the triple-gluon correlation functions 
based on the preliminary data on A® by RHIC. 



1. Introduction 

The single spin asymmetry (SSA) in inclusive hard processes appears as a twist-3 observable 
in the collinear factorization which is valid in describing the large-Pp hadron productions. In 
this twist-3 mechanism, SSA is represented in terms of multi-parton correlations in the hadrons. 
Such correlations in the transversely polarized nucleon can be represented by the quark-gluon 
correlation functions and the triple-gluon correlation functions. So far the former effect has 
been investigated for various processes. To probe the latter effect, open-charm production in 
ep and pp collisions is an ideal tool, since the cc pair is created mainly by the gluon-photon or 
gluon-gluon fusion processes. 

In our recent paper [I] , we formulated a method of calculating the contribution of the triple- 
gluon correlation functions to the single-spin-dependent cross section for the D-meson production 
in semi-inclusive deep inelastic scattering (SIDIS), ep 1 ' — > eDX. There we identified the complete 
set of the relevant triple-gluon correlation functions and derived the corresponding cross section 
in the leading order with respect to the QCD coupling constant. Our result differed from the 
previous study in the literature [2] , and we have clarified the origin of the discrepancy. 

In this work, we will apply the formalism of pQ to the pp-collision, p^p — > DX, and derive the 
corresponding single-spin-dependent cross section. This study is relevant to the ongoing RHIC 
experiment. We also provide a simple estimate for the constraint on the triple-gluon correlation 
functions, using a preliminary data on p^p — > DX reported at RHIC [3]. 



2. Triple-gluon correlation functions for the transversely polarized nucleon 

Triple-gluon correlation functions in the transversely polarized nucleon are denned as the 
correlation functions of the three gluon's field strength tensor F a ^ . Since there are two ways 
of constructing color-singlet correlation functions from the contraction with the symmetric and 
antisymmetric structure constants for the color SU(3) group, d\, ca and fb ca , one can define two 
gauge-invariant correlation functions in the nucleon as [T] 
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where Mn is the nucleon mass, S is the transverse-spin vector for the nucleon, n is the light-like 
vector satisfying p ■ n = 1 and we used the shorthand notation as F^ n = F /3p n p etc. Hermiticity, 
PT-invariance and permutation symmetry lead to the decomposition of ([1]) and ([2]) in terms of 
the two real functions 0(x\,x 2 ) and N(x\,x 2 ) which have the following symmetry properties, 

0(xi,x 2 ) = 0(x 2 ,xt), 0(xi,x 2 ) = 0(-Xi,-x 2 ), 
N(x 1 ,x 2 )=N(x 2 ,x 1 ), N( Xl ,x 2 ) = -N(- Xl ,-x 2 ). (3) 

The gauge-link operator which restores gauge invariance of the correlation functions is suppressed 
in ([!]) and ([2]) for simplicity. 
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3. Polarized cross section formula for p^p — > DX 

The formalism for calculating the contribution of the triple-gluon correlation functions to 
ep^ —> eDX developed in pQ can be directly applied to p*p — > DX [4j. In this process, a cc-pair 
is created by the gluon-fusion process, and the initial-state-interaction (ISI) diagrams as well as 
the final-state-interaction (FSI) diagrams give rise to the single-spin-dependent cross section as 



a pole contribution at x\ = x 2 . The twist-3 cross section for p^(p, S±) + p(p') 
with the center-of-mass energy yfs is given by 
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where 5 C = 1 and <5g = —1, Df(z) represents the c — > D or c — > D fragmentation functions, G(x') 
is the unpolarized gluon density, and p c is the four-momentum of the c (or c) quark (mass m c ) 
fragmenting into the final D (or D) meson. The partonic hard cross sections in are given by 

( l 2 + u 2 + 4m 2 s 4mis2 
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where N = 3 and C F = (N 2 - 1)/(2N), and s, t, u are denned as 

s = (xp + x'p') 2 , t = (xp — p c ) 2 — m 2 , u = (x'p' — p c ) 2 — m 2 . (7) 

The hard cross sections associated with the first term in the first parentheses in (0) and ^ come 
from ISI, and those associated with the second and the third terms in the same parentheses come 
from FSI. As in the case of ep^ — > eDX, the cross section in (jU) receives the contribution from 
the four functions 0(x,x), O(x,0), N(x,x), N(x,0). Unlike the case of SIDIS, presence of ISI 
gives rise to the different hard cross sections for O and N functions. From (j3J), it is clear that the 
process p^p — > DX itself is not sufficient for the complete separation of the four functions. For 
the separation, the process ep^ — > eDX serves greatly, since it has five structure functions 
with different dependences on the azimuthal angles to which the four functions contribute 
differently pQ . 

Our result in @ differs from a previous work [5]: The result in [5] is obtained from @ 
by omitting the terms with 0(x, 0) and N(x, 0) (and their derivatives) and by the replacement 
0(x, x) — > 0(x, x) + 0(x, 0) and N(x, x) — > N(x, x) — N(x, 0). This difference originates from an 
ad- hoc assumption in the factorization formula in [5l[2]. We emphasize the appearance of the four 
different contributions with {0(x, x), 0(x, 0),N(x, x), N(x, 0)} is a consequence of the symmetry 
property implied in the decomposition ([1]) and ([2]), in particular, the different coefficient tensors 
in front of 0(x, x) and 0(x, 0) (likewise for N(x, x) and N(x, 0)) at x\ = X2 = x lead to different 
hard cross sections for the above four functions. See [H 0] for more details. 

4. Numerical estimate 

As is shown in @, four nonperturbative functions are involved in the twist-3 cross section for 
^4^. At present there is no information on these functions. Preliminary data on A® by RHIC- 
PHENIX [3] suggests \Ajy\ < 5 %. So we will present a model estimate on the upper bound of 
these nonperturbative functions for several cases by requiring that the calculated A® be less 
than 5 %. We first found from © and © the relations l^ 01 ' 02 ^ 1 ^ 2 ! > |^03,04,iV4,iV4| and 
a 01 ~ g° 2 ~ a N1 ~ —a N2 for the RHIC kinematics (see below). Accordingly, if 0(x,x) and 
O(x,0) have the same (opposite) sign, they contribute to A® constructively (destructively). 
Likewise for N(x,x) and A^(x,0). From this fact we present a model calculation of A® for two 
extreme cases: 

Model 1 : 0(x, x) = -0(x, 0) = N(x, x) = N(x, 0), (8) 
Model 2 : 0(x, x) = 0(x, 0) = N(x, x) = -N(x, 0). (9) 

In Model 1 the contributions to A® from four functions become minimum, while in Model 2 
they become maximum. Therefore larger magnitude for the functions is allowed for Model 1 in 
order to make |^4^| < 5 %. For the actual form of the nonperturbative functions, we follow [2 
and set 
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where G(x) is the unpolarized gluon distribution and Kq is a constant which we determine 
so as to achieve \A N \ < 5 % for the RHIC kinematics. Obviously, the ansatz (|10p is a very 
crude approximation and the result below should be taken only as an estimate of the order-of- 
magnitude. For the numerical calculation, we use GJR08 distribution [6j for G(x) and KKKS08 
fragmentation function [7] for Df{z). We also assumed the same scale dependence for 0(x,x) 
etc as G(x) for simplicity. We calculated An for the D and D mesons at the RHIC energy of 
yfs = 200 GeV and Pt = 2 GeV with the parameter m c = 1.3 GeV by setting the scale of all 

the distribution and fragmentation functions at \i = \JPj< + ni*. 

Fig. 1 shows the result of A® for Model 1 with Kq = 0.005. In order to see the behavior of 
each term in Q , we showed in the left figure of Fig. 1 the contribution to A N ° from the four 
functions 0(x,x), O(x,0), N(x,x) and N(x, 0). (For N(x,0), —A N is shown.) As is seen from 
this figure, if one sets all four functions identical, contributions to A N ° from them are very close. 
In the middle and right figures of Fig. 1, A N for the D° and D° mesons are shown, respectively. 
Even though each of four contribution gives rise to as large as 50 % asymmetry as shown in the 
left figure, the resulting total An is of 0(5%) due to the cancellation among them. Therefore 
if the relation ([8]) approximately holds, Kq = 0.005 in (fTQ|) provides an upper limit for 0(x,x) 
etc. Because of the relation ((8|), the O-term and the N-term contribute to An constructively 
(destructively) for D° (D°) meson as shown in Fig. 1. This feature was also observed in [5]. 




Xp Xp Xp 



Figure 1. A N for Model 1. Decomposition of A N into the 4 components (left). An for D' 
meson (middle). An for D° meson (right). 



Fig. 2 shows the result for Model 2 with Kq = 0.0001. In the left figure of Fig. 2, the 
behavior of each contribution is shown. In the middle and right figures, A N for the D° and D° 
mesons are shown, respectively. In this Model 2, we had to take Kq as small as Kq ~ 0.0001, 
since 0{x,x) and O(x,0) as well as N(x,x) and iV(x,0) contribute to An of either D° or D° 
constructively. 

From the left figures of Figs. 1 and 2, one can see that the magnitude of each of the four 
contributions is very similar if the four functions have the same magnitude. This means that 
the contributions to A N from 0(x,x) + O(x,0) and N{x,x) — N(x,0) are, respectively, much 
larger than those from 0(x, x) — 0{x, 0) and N{x, x) + N(x, 0), because 

\<y + o |3>|cr —a |, |cr —a |3>|<r + a \. (11) 

Since there is no reliable nonperturbative information on these functions, it is natural to expect 
that all four functions 0(x, x) ± 0(x, 0) and N(x, x) =F N(x, 0) have a similar magnitude. In this 
case, A N is mostly determined by 0(x, x) + 0(x, 0) and N(x, x) — N(x, 0), while 0(x, x) — 0(x, 0) 
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Figure 2. A N for Model 2. Decomposition of A N into the 4 components (left). An for D 
meson (middle). An for D° meson (right). 



and N(x, x) + N(x, 0) can be neglected by the relation (jlip . From this observation, one obtains 
a modest estimate on the upper bound for the combinations as 

\0{x, x) + 0(x, 0)| < (0.0003 ~ 0.0004)xG(x), 

\N(x, x) - N(x, 0)| < (0.0003 ~ 0.0004)a;G(x), (12) 

if \A N \ < 0(5%). We remind that the relation (jlip is a peculiar feature for the process 
p^p — > DX and does not generally hold in ep^ — > eDX pQ . 

5. Summary 

In this work, we have studied the contribution of the triple-gluon correlation functions to A N 
in p^p — > DX. We derived the corresponding twist-3 single-spin-dependent cross section in the 
leading order with respect to the QCD coupling constant. The complete cross section receives 
contribution from the four functions 0(x, x), O(x,0), N(x,x) and N(x,0) as in the case of 
ep^ — > eDX pQ , and differs from the result in a previous work [5] . We have also presented a 
model calculation for A N at the RHIC energy. Assuming \ A N \ < 0(5%) as suggested by the 
RHIC preliminary data, we obtained a modest estimate (|12p for a particular combination of the 
triple-gluon correlation functions. The detail of our analysis will be reported elsewhere 
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